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Abstract The NEMO Collaboration is conducting an 
R&D activity towards the construction of a Mediter- 
ranean km 3 neutrino telescope. In this work, we present 
the results of Monte Carlo simulation studies on the ca- 
pability of the proposed NEMO telescope to detect and 
identify point-like sources of high energy muon neutri- 
nos. 
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1 Introduction 

The detection of high energy neutrinos is considered one 
of the most promising means to investigate non-thermal 
processes in the Universe. A first generation of small 
scale detectors has been realized (NT-200 pQ in the Baikal 
lake and AMANDA at the South Pole), demonstrat- 
ing the possibility to use the Cerenkov technique to track 
high energy neutrinos. Besides, these experiments have 
set limits on neutrino fluxes. Other small scale detectors 
are at different stage of realization (ANTARES [3j and 
NESTOR |U). 

Actual expectations on neutrino fluxes, mainly based 
on the measured cosmic ray fluxes and the estimated 
fluxes for several high energy sources from theoretical 
models 0, require detectors of km 3 scale. Following the 
success of AMANDA, the largest operating detector, the 
realization of the IceCube km 3 detector j^j has started 
at the South Pole. On the other hand, many issues, as 
the full sky coverage, strongly support the construction 
of a km 3 scale detector in Mediterranean Sea. 
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2 The NEMO project 

The NEMO Collaboration J7J is performing R&D to- 
wards the design and construction of the Mediterranean 
km 3 neutrino detector. The activity was mainly focused 
on the search and characterization of an optimal site 
for the detector installation and on the development of a 
feasibility study for the detector. A deep sea site with op- 
timal features in terms of depth and water optical prop- 
erties has been identified at a depth of 3500 m about 
80 km off-shore Capo Passero and a long term monitor- 
ing of the site has been carried out 0. The feasibility 
study of the km 3 detector includes the analysis of all the 
construction and installation issues and the optimization 
of the detector geometry by means of numerical simula- 
tions. The validation of the proposed technologies via an 
advanced R&D activity, the prototyping of the proposed 
technical solutions and their relative validation in deep 
sea environment will be carried out with the two pilot 
projects NEMO Phase- 1 and Phase-2 0. 



3 Detector lay-out 

The geometry of the NEMO-km 3 telescope simulated in 
this work is a square array of 9 x 9 towers equipped with 
5832 optical modules (10" diameter PMTs) 0. The tow- 
ers, moored on the seabed at 3500 m depth, have an 
instrumented height of 680 m, with a storey-storey dis- 
tance of 40 m, and are spaced of 140 m. The detector 
response is simulated using the codes developed by the 
ANTARES Collaboration [THj. In the simulation codes, 
the light absorption length, measured in the site of Capo 
Passero (L a w 68 m at 440 nm [H]), is taken into ac- 
count. Once the sample of PMT hits is generated, spuri- 
ous PMT hits, due to the underwater optical noise ( 40 K 
decay), are introduced, with a rate of 30 kHz for 10" 
PMTs, corresponding to the average value measured in 
Capo Passero site. 

The simulated detector lay-out reaches an effective 
area of 1 km 2 at a muon energy of about 10 TeV and an 
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angular resolution of a few tenths of degrees at the same 
energy as plotted in Fig. 



gular resolution a = 0.19° ± 0.02° (see Fig.EJ. Besides 
detecting its position in the sky allows us to determine 
the absolute orientation of the detector • 



(J) 
U 

cr 
c 

1 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 
0, 

1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




' I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ' 



2.5 3 3.5 



4.5 5 5.5 6 

Log(s M / GeV) 



■ ■ ■ ■ ' ■ ■ ■ ■ i ■ ■ i ■■■ ■ ' ■ ■ ■ ■ I ■ ■ ■ ■ ' ■ ■ ■ ■ 



J., 



2.5 3 3.5 4 4.5 5 5.5 6 

Log(c M / GeV) 

Fig. 1 Effective area and angular resolution of the simulated 
NEMO detector as a function of the muon energy. The detec- 
tor response is simulated considering a diffuse flux of up-going 
muons reaching the detector surface. 



4 Detector pointing accuracy 

The detector angular resolution is one of the most impor- 
tant parameters in the identification of point-like sources. 
Therefore, an experimental determination is required. A 
possible method, already adopted in cosmic ray detec- 
tors, consists in the observation of the so called Moon 
shadow Since the Moon absorbs cosmic rays, we ex- 
pect a lack of atmospheric muons from the direction of 
the Moon disk. The detection of the muon deficit pro- 
vides a measurement of the detector angular resolution. 

Monte Carlo simulations show that the NEMO tele- 
scope could be able to detect the Moon shadow and that 
about 100 days are needed to observe a 3a effect. As- 
suming a detector point spread function with a Gaus- 
sian shape, the same simulations yields a detector an- 
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Fig. 2 Detected muon event density versus the angular dis- 
tance from the Moon center, assuming 1 year of data tak- 
ing. The points are fitted assuming a Gaussian shape for 
the point spread function, obtaining an angular resolution 
of (7 = 0.19° ±0.02° H2. 



5 NEMO-km 3 sensitivity to neutrinos from 
point-like sources 

5.1 Calculation of the detector sensitivity 

The detector sensitivity spectrum is calculated according 
to the following formula: 
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where ~pg (b) is the 90% c.l. average upper limit for an 
expected background b and calculated as suggested by 
Feldman and Cousins |13j: (depu/ds^o is an arbitrary 
point source spectrum inducing a mean signal N^. The 
detector sensitivity is calculated taking into account both 
atmospheric neutrino and muon backgrounds. 



5.2 Simulation of the atmospheric muon and neutrino 
background 

A sample of 7- 10 9 atmospheric neutrinos have been gen- 
erated using the ANTARES event generation code, based 
on a weighted generation technique |10j . The events were 
generated in the energy range 10 2 -j- 10 8 GeV, with a 
spectral index X = 2 and a 4-7T isotropic angular distri- 
bution. The events were then weighted to to the sum of 
the Bartol flux and of prompt neutrino rqpm model 
[T3] flux. When the event weight is calculated, the neu- 
trino absorption in the Earth, as a function of neutrino 
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energy and direction, is taken into account jXG . So doing, 
we compute a number « 4 • 10 4 of detected atmospheric 
neutrino events per year of data acquisition. 

Atmospheric muons are generated at the detector, 
applying a weighted generation technique. We generated 
a sample of N to tai — 2.5 • 10 7 muons, in the energy range 
1 TeV -7- 1 PeV, with a generation spectral index X = 3. 
We also generated N tota i = 4- 10 7 events in the range 100 
GeV - 1 TeV, with a generation spectral index X = 1. 
Muons are generated with an isotropic angular distribu- 
tion. The events are weighted to the Okada parameter- 
ization |TJ|, taking into account the depth of the Capo 
Passero site (D = 3500 m) and the flux variation inside 
the detector sensitive height (h « 900 m). According 
to the Okada parameterization, the expected number of 
reconstructed muon events is about 4 • 10 8 per year. At- 
mospheric muon parameterization by Bugaev et al |18j 
is also considered. In this case the expected number of 
reconstructed events is 5 • 10 s per year but no signifi- 
cant differences are observed in the detector sensitivity 
values. 

The simulated statistics cover only a few days. Con- 
sidering that reconstructed events have a flat distribu- 
tion in Right Ascension (RA), we can project the simu- 
lated events in a few degrees bin Z\RA, centered in the 
source position. So doing, we get statistics of atmospheric 
muons corresponding to a time ^1 year at all source de- 
clinations. 

Major details on the Monte Carlo simulation of the 
NEMO telescope could be found in ref. |lfi| . 

5.3 Criteria for the atmospheric background rejection 



5.4 Detector sensitivity to neutrino point-like sources 

In this section, we calculate the expected detector sensi- 
tivity to neutrinos from point-like sources. We simulated 
muons induced by ~ 10 9 neutrinos with energy range 
10 2 -hl0 8 GeV and X=l. These events are weighted to the 
neutrino spectrum ((kp v /de v )a = 10' 7 £~Q e y (GeV -1 
cm' 2 s' 1 ). As a first case, we consider a source at a 
declination of 8 = —60°. Such a source has a 24 hours of 
diurnal visibility and it covers a large up-going angular 
range (d^ = 24° — 84°). The sensitivity for this source 
is therefore representative of the average response of the 
NEMO detector. 

In Table Q we report the expected sensitivity for 
different values of the spectral index a, considering 3 
years of data taking. In Fig. |3| we show how the en- 
ergy spectrum of reconstructed events is reduced after 
the event selection, varying the neutrino spectral index 
a. The same plots show also as the spectrum peak moves 
towards lower energies for softer spectral indices. 



Table 1 Sensitivity to a point-like neutrino source at S — 
—60°, for different spectral indices a and 3 years of data 
taking. The sensitivity spectrum e"(dtp v ld£ v )$o is expressed 
in GeV Q ~V cm 2 s. 
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The used reconstruction algorithm is a robust track fit- 
ting procedure based on a maximization likelihood me- 
thod 10 . In this work, we used, as a goodness of fit 
criterion, the variable: 

yl __tog Y C) +0 . 1(iVcomj) _ 1 ) ) (2 ) 
Ndof 

where log(£) /TVdof is the log- likelihood per degree of 
freedom (Ndof) and N comp is the total number of com- 
patible solutions found by the reconstruction program. In 
particular, events are selected if the variable A is greater 
than a given value A cut . This quality cut is here applied 
together with other selection criteria as listed in the fol- 
lowing: 

— the number of hits Nja , used to reconstruct the muon 
track, must be greater than a given value Nf$\ 

— the muon must be reconstructed with ?9'' ec < $™ ax 
in order to reject down-going events; 

— only events reconstructed in a circular sky region cen- 
tered in the source position and having a radius of 
rfcn are considered. 

The optimal values of A cut , Nj^, $™ ax and r^n are 
chosen to optimize the detector sensitivity. 



The detector sensitivity was calculated as a function 
of the years of data taking and reported in Table [21 Re- 
sults for a — 2 are also plotted in Fig. 01 compared to the 
IceCube sensitivity obtained for a 1° search bin [Hj. Our 
results show that the proposed NEMO detector reaches a 
better sensitivity to muon neutrino fluxes with a smaller 
search bin. 



Table 2 Sensitivity to a point-like neutrino source at 5 = 
—60°, for different spectral indices a and for different 
number of years of data taking. The sensitivity spectrum 
eZ{dip v / dE v )gQ is expressed in GeV Q_1 / cm 2 s. 
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Fig. 3 Neutrino energy spectra of reconstructed events, be- 
fore (solid line) and after (dashed line) the event selection. 
The spectra are normalized to the sensitivity fluxes for 3 years 
of data taking. 
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Fig. 4 Sensitivity to a neutrino spectrum with a = 2, coming 
from a 8 — —60° declination point-like source and comparison 
with the IceCube detector El. 



ity doesn't get worse until that e™ ax reaches the energy 
peak of reconstructed neutrino spectra (see Fig. |3J. In 
the case of hard spectrum sources, the detector sensitiv- 
ity is better and it gets better if the spectrum extends 
to VHE. Expected astrophysical neutrino spectra could 
be softer (a^2); in this case the sensitivity value doesn't 
vary much with e™ ax . 

We finally consider the dependence of the sensitiv- 
ity on the source declination. In particular, the detector 
sensitivity gets worse with increasing declination due to 
the decrease of the diurnal visibility (to respect with the 
latitude of the Capo Passero site) . Fig. shows the sen- 
sitivity versus the source declination, considering three 
years of data taking and a — 2. The worst sensitivity is 
2.5 • 10~ 9 GeV/cm 2 s, calculated for a source declination 
of S = 50° for which the diurnal visibility reduces to a 
few hours per day. 



6 Physics cases 

In this section, we consider the case of two particular 
sources: microquasar LS 5039 and SNR RX J1713.7- 
3946, both observed in the TeV gamma-ray region. For 
each source, we compute the detector sensitivity and the 
expected number of source events compared with the 
background. 



6.1 Microquasar LS 5039 



The expected astrophysical neutrino spectra could 
not extend up to 10 8 GeV, especially in the case of 
Galactic sources. For this reason, we also computed the 
detector sensitivity as a function of the high energy neu- 
trino cut-off e" lax . Results of our calculations are plot- 
ted in Fig. Decreasing the energy cut-off, the sensitiv- 



The H.E.S.S. telescope has recently detected TeV 7-rays 
from LS 5039 This discovery provided the first un- 
ambiguous evidence for presence of multi-TeV particles 
in microquasars. 

Aharonian et al. [20] discussed different possible sce- 
narios for the production of the observed 7-ray flux. They 



Detection potential to point-like neutrino sources with the NEMO-km 3 telescope 



5 



3.5 



10 



CO 



E 10 
o 

-No" 7 

a 

> -8 
45 10 

o 

o ,,-r 9 
^ 10 

to 

^10 t- 

"O 10 



co -12 
10 



-13 

10 b- 

-14 

10 t- 



10 



-15 





a=2.5 


1 1 


— ¥ — ▼ — ▼ — ▼ — ▼ 




a=2 


1 1 IIMIIj 

1 


— * — * — a — * — ▲ 




a=1.5 
^* ■ — ■— ■— ■ 


1 1 IIMIIj 1 








1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 iTi 1 1 1 



.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 

"™ 7 GeV) 



Log(£„ 



£ 
o 

> 

o 



3 - 



2.5 



S 1-5 



S 
^3 



0.5 



_i i i I i i i_ 



-80 



-60 



-40 



-20 



20 



. I . . 

40 60 

<5 (deg) 



Fig. 6 Sensitivity to a point-like neutrino source at declina- 
tion S, for spectral index a — 2 and 3 years of data taking. 

Table 3 Detector sensitivity to neutrinos from LS 5039. The 
sensitivity /^,9o is expressed in units of erg/cm 2 s and refers 
to a detector live time of 3 years. The corresponding values 
of i!}™ ax , A C ut and run are also given. We also report the 
number N™ from the source compared to the atmospheric 

background events N*, surviving the event selection. 
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Fig. 5 Sensitivity to a point-like neutrino source at 8 = 
—60°, for different spectral indices a and for different val- 
ues of the high energy cut-off. 



considered both leptonic and hadronic production mech- 
anisms and argued in favor of a TeV photon flux orig- 
inating from pp interaction. If so, 7-rays should be ac- 
companied by TeV neutrinos with an average energy flux 
of fj 1 = 1(T 10 erg/cm 2 s. 

In Tab. are summarized the detector sensitivities 
for microquasar LS 5039, assuming neutrino fluxes with 
spectral indices a = 1.5 and 2, in the energy range 0.1 
TeV and e™ ax = 10 and 100 TeV. The expected number 
of selected events, induced by the flux is given in 
the same table. The comparison with the atmospheric 
background shows that an evidence could be expected in 
a few years of data taking. 



6.2 SNR RX J1713. 7-3946 

The CANGAROO Collaboration observed 7-rays from 
SNR RX J1713. 7-3946, claiming the hadronic origin of 
the measured energy spectrum |21| . Alvarez-Muniz and 
Halzen (A&H) [221 calculated the high -energy neutrino 



flux associated with this source. Their calculations yield 
to an expected neutrino spectrum 



eZ{dLp v /de u ) th = 4.14 • 10~ 8 cm- 2 s" 



GeV" 



(3) 



with spectral index a = 2 and extending up to ~ 10 TeV. 
More recent calculations were performed by Costantini 
and Vissani (C&V) [22] > based on the 7-rays flux mea- 
sured by the H.E.S.S. experiment |24|. According to their 
calculations, we expect a neutrino spectrum 
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with a — 2.2 and a neutrino energy ranging between 50 
GeV and 1 PeV. 

In this paper, we calculate the expect detector sen- 
sitivity for RX J1713. 7-3946. Despite it is an extended 
source, in a first approximation we can consider it as 
point-like since its diameter (0 = 1.3°) [23 is compara- 
ble o smaller than the detector search bin. Results are 
reported in Tab. 01 where the dependence of the sensitiv- 
ity on the spectral index a and on the high energy cut-off 
£ max - 1S s i 10wn> j n the- same table, we also report the num- 
ber of expected events surviving the selection criteria, 
considering both the theoretical predictions. The com- 
parison with the atmospheric background shows that, 
also in this case, the NEMO telescope could identify the 
source in a few years of data taking. 
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Table 4 Detector sensitivity to neutrinos from the SNR RX 
J1713. 7-3946. The sensitivity is expressed in units of cm -2 
s _1 GeV -1 and refers to a detector live time of 3 years. The 
corresponding values of $™ ax , A cu t and Thin fire also given. 
We also report the number N™ from the source compared to 

the atmospheric background events surviving the event 
selection. 
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7. The NEMO web site: nemoweb.lns.infn.it. 
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7 Conclusions 

The possibility to detect TeV muon neutrinos from point- 
like sources with the proposed NEMO-km 3 underwater 
Cerenkov neutrino telescope has been investigated. In 
particular Monte Carlo simulations were carried out to 
determine the expected response of the km 3 telescope. 

Our simulations show that it could be possible to 
observe the Moon shadow in about 100 days. This de- 
tection provides a measurement of the detector angular 
resolution and of the detector pointing accuracy. 

We also computed the detector sensitivity to muon 
neutrinos from point-like sources, defined as the min- 
imum flux detectable with respect to the atmospheric 
muon and neutrino background. The dependence of the 
sensitivity on the neutrino spectral index and energy 
range, on the source declination and on years of data 
taking has been studied. 

Finally, we consider the case of two particular sources: 
microquasar LS 5039 and SNR RX J1713. 7-3946; both 
observed in the TeV gamma-ray region. For each source, 
we compute the detector sensitivity and the expected 
number of source events compared with the background. 
Our results show that, assuming present predictions of 
TeV neutrino fluxes, the proposed NEMO telescope could 
identify both sources in a few years of data taking. 
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